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ABSTRACT 


A  key  to  modeling  polarimetric  signatures  of  painted  objects  is  the  polarimetric  bidirectional  reflectance 
distribution  function  (BRDF).  In  particular,  it  is  important  to  correctly  capture  the  out-of-plane  behavior 
of  the  BRDF.  From  the  theoretical  point  of  view,  the  BRDF  is  a  complex  object  -  a  4  x  4  matrix  of 
functions,  each  depending  on  three  angles.  Brute  force  tabulation  of  such  an  object  is  unlikely  to  ever  be 
a  feasible  modeling  approach.  For  this  reason,  model  forms  of  the  BRDF,  rooted  in  physics,  are  very 
valuable  for  extending  measured  data  and  for  use  in  modeling  codes.  Only  a  limited  number  of  such 
models  have  been  mentioned  in  the  literature.  Perhaps  the  one  most  relevant  to  paints,  such  as  CARC 
(Chemical  Agent  Resistive  Coating)  coatings,  is  the  microfacet  model.  Unpolarized  versions  of  this 
model  have  been  used  in  various  forms  for  years.  This  paper  presents  a  full  polarimetric  version  of  this 
model  and  discusses  its  properties.  Recently  acquired  polarimetric  data  is  compared  with  the  model.  It  is 
found  that  the  model  is  successful  in  capturing  qualitative  aspects  of  this  data. 


1.  Introduction 


An  important  issue  in  polarimetric  target  signature  modeling  is  the  form  of  the  polarimetric  bidirectional 
reflectance  distribution  function  (BRDF)'  ’^.  The  general  form  of  the  polarimetric  BRDF,  a  4  x  4  matrix 
function,  is  very  complex  with  subtle  constraints  linking  the  sixteen  elements.  Since  quantities  such  as 
the  degree  of  polarization  involve  the  difference  of  large  numbers  (to  form  a  small  number),  details  of  the 
form  of  the  BRDF  are  potentially  quite  important.  This  is  especially  true  for  out-of-plane  scattering 
geometries.  There  is  a  clear  advantage  in  working  with  a  closed  form  expression  which  derives  from  a 
physical  model,  exhibits  reciprocity  symmetry,  and  integrates  correctly  to  the  directional  reflectivity. 
Calculations  using  such  a  BRDF  will  be  free  of  artifacts  that  could  arise,  for  example,  from  using  an 
interpolation  of  laboratory  measurements  not  obeying  the  constraints  among  the  sixteen  elements.  This 
paper  presents  the  closed  form  result  for  the  polarimetric  BRDF  derived  from  the  microfacet  model. 
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The  microfacet  model  is  rooted  in  geometric  optics.  Various  forms  of  this  model  have  been  in  the 
literature  since  the  1960s.^'^  In  its  simplest  form  the  model  envisions  a  collection  of  randomly  oriented 
microfacets  as  comprising  the  rough  surface.  Each  microfacet  acts  as  a  specular  reflector.  A  statistical 
distribution  of  facet  orientations  described  by  a  distribution  function  is  postulated.  Only  a  reflection  law 
for  the  microfacets,  such  as  that  provided  by  the  Fresnel  relations,  is  needed  to  complete  a  definition  of 
the  model.  The  literature  contains  several  derivations  of  the  scalar  (unpolarized)  BRDF  arising  from  this 
model.  These  treatments  use  a  variety  of  approaches  and  notations.  Rather  than  review  this  derivation 
here,  the  scalar  result  is  simply  quoted. 

An  important  relation  links  the  unpolarized  BRDF,/  to  the  directional  hemispherical  reflectivity  (DHR), 

DHR{e,)=  (i) 

The  BRDF  is  a  function  of  the  incident  polar  angle  9\,  the  scattered  polar  angle  0^,  and  the  difference 
between  the  incident  and  scattered  azimuthal  angles,  where  we  assume  the  sample  is  rotationally 
symmetric.  The  solid  angular  integration  is  over  the  hemisphere  above  the  surface.  The  integration 
differential,  6^ (O, ,  is  given  by  sin(0r)  d^rd^r .  Further,  the  scalar  BRDF  must  be  an  even  function  of 
The  cosine  in  the  integrand  appears  because  the  BRDF  is  defined  as  being  per  unit  projected  area.  Eq.  (1) 
is  a  normalization  condition  on  the  BRDF.  An  ideal  Fambertian  reflecting  surface  has  the  property  that 
the  outgoing  power  per  unit  solid  angle,  per  unit  projected  area,  is  independent  of  angle.  This  leads  to  the 
result  that /=  l/n,  that  is,  a  constant,  for  the  Fambertian  case.  Another  important  result  concerns  the 
radiance  reflected  from  a  horizontal  surface  exposed  to  an  incoming  sky  radiance  profile  Rsky(^sky)-  The 
definition  of  the  BRDF  leads  to  the  result 

^reflected  r  )  “  ^skyfi^sky  ^sky  )^sky  sky  )  sky  )  '  (2) 

This  result  generalizes  immediately  to  the  case  where  the  sky  radiance  depends  on  the  azimuthal 
coordinate.  It  is  the  main  result  needed  for  modeling  reflected  signatures.  Consideration  of  the  process 
reciprocal  to  that  leading  to  the  expression  for  the  DHR  leads  to  a  reciprocity  condition  on  the  BRDF. 
These  can  be  written  as 


\mAA -^,)cos{e^)d^co^  =  RHO{e^) 


(3) 


where  d^<y  =  sin(^dfti^.  The  functional  forms  on  the  right-hand  sides  of  Eqs.  (3)  are  the  same  for  both 
integrals  and  is  the  functional  form  of  the  DHR  .  A  sufficient  condition  on / to  insure  this  reciprocity 
condition  is 


(4) 

As  will  be  seen  below,  the  microfacet  model  assumption  leads  to  a  form  obeying  this  interchange 
symmetry. 

The  BRDF  is  closely  related  to  another  object,  a,  which  appears  in  some  of  the  older  literature^.  This 
object  is  the  differential  scattering  cross  section  per  unit  illuminated  area.  The  definition  of  a  differential 
scattering  cross  section  implies  that  the  power  per  unit  solid  angle  going  into  the  viewing  direction  is  the 
product  of  the  incident  flux  and  the  differential  scattering  cross  section.  It  can  be  shown  that  the  BRDF 
is  related  to  <7  by  a  simple  relation. 


a 


(5) 


/  = 


cos(0Jcos(0,.) 


2.  Definitions  of  angies 


Some  angular  definitions  were  introduced  in  Section  1.  This  section  collects  all  the  notation  used  to  refer 
to  angles  appearing  in  the  sections  below. 

Three  sets  of  spherical  polar  angles  and  associated  unit  vectors  will  prove  to  be  useful.  These  are:  1)  a 
unit  vector  ri ,  and  associated  polar  and  azimuthal  angles  9\  and  pointing  to  the  direction  of  the  source 
radiation;  2)  a  unit  vector  Fr ,  and  associated  angles  6,  and  (pr,  pointing  toward  the  position  of  a  viewer;  3) 
a  unit  vector  n,  and  associated  angles  9  and  (j),  pointing  in  the  direction  of  the  local  surface  normal  of  a 
microfacet.  For  each  of  these  three,  the  symbol  co  (with  appropriate  subscript  in  the  case  of  the  first  two) 
will  be  used  to  refer  to  the  polar,  azimuth  pair.  The  corresponding  two-dimensional  differential  element 
is  d^<y  =  sin(^  d^d^.  The  overall  surface  normal  will  be  denoted  as  z  and  this  direction  will  serve  as  the 
polar  axis  for  the  spherical  polar  coordinates.  A  useful  auxiliary  angle  is  P,  the  angle  of  incidence  onto 
or  angle  of  reflection  from  a  microfacet.  In  the  microfacet  model  these  are  the  same  angle,  and  we  have 

cos{P)  =  n  ■  r,.  =  n  ■  .  (6) 

Also  following  from  the  assumed  mirror- like  nature  of  a  microfacet  is  the  vector  relation 

r,.  +  =  2  cos(/?)n  .  (7) 

Useful  equations  can  be  obtained  from  this  vector  relation  by  forming  dot  products  with  both  sides.  For 
example,  taking  the  dot  product  of  Eq.  (7)  with  z  leads  to 


cos(0)  = 


cos(0,)-l-cos(0^) 

2cos(y5) 


3.  Formal  definition  of  the  microfacet  model 


(8) 


The  microfacet  model  assumes  that  a  collection  of  small  microfacets  comprises  the  rough  surface.  Each 
microfacet  is  a  specular  reflector  obeying  Snell’s  law  of  reflection  (angle  of  incidence  equal  angle  of 
reflection)  with  reflectivity  given  by  a  Fresnel  reflectivity’,  p{P),  depending  on  the  local  angle  of 
incidence,  p.  Each  microfacet  is  characterized  by  a  normal  unit  vector  n,  and  the  microfacets  are  assumed 
to  be  distributed  symmetrically  about  z  in  orientation  according  to  a  distribution.  A  well-studied  example 
is  the  Gaussian  form. 


p{e)  = 


1 


2na  cos  (0) 


-tan’(6>)^ 


(9) 


This  is  a  normalized  form,  so  that 
jd^(Up(9)  =  l, 


(10) 


where  the  integral  is  over  the  upper  hemisphere.  The  local  surface  slope  is  tan(^.  The  slope  variance  is 
o^.  This  form  can  be  shown  to  arise  from  the  assumption  of  a  Gaussian  surface  height  distribution,*  and 
has  the  advantage  of  enforcing  0<  7i/2. 


4.  The  scalar  BRDF 


Here  the  result  for  the  scalar  BRDF  is  quoted  without  derivation.  The  scalar  BRDF  is 


,3-6 


f- 


1  1  1  2(j" 


exp(- 


tan"(6>) 


) 


2n  4(j^  cos"^(0)  cos(0Jcos(0.) 

In  working  with  Eq.  (1 1),  we  can  use: 

cos(0,.)  +  cos(0J 


piP)- 


(11) 


COS(0)  : 


2cos(/?) 


(12) 


cos(2y5)  =  cos(0, ) cos(0^ )  +  sin(0. ) sin(0^ ) cos(^,. 


with  the  former  being  a  restatement  of  Eq.  (8)  and  the  latter  expression  coming  from  the  spherical 
addition  theorem. 

With  these  auxiliary  formulas,  it  can  be  seen  that  the  above  expression  for /has  the  interchange  symmetry 
discussed  in  Section  1  and  thus  obeys  the  reciprocity  condition.  This  expression  agrees,  apart  from  a 
factor  of  471,  with  the  result  for  the  differential  scattering  cross  section  derived  in  Ref  3.  The  factor  of  471 
results  from  the  way  the  differential  scattering  cross  section  is  defined  and  normalized  in  Ref  3. 

5.  Properties  of  the  BRDF 


In  addition  to  the  interchange  symmetry,  we  note  here  another  property  of  the  microfacet  BRDF.  If 
piji)  =  1 ,  then 

=  DHR{0,)  =  1  (13) 


for  any  value  of  0.  The  integration  over  (O^  is  rather  involved  because  of  the  non-trivial  out-of-plane 
behavior  of f.  Rather  than  attempt  this  integration  analytically,  the  commercial  program  MATHCAD 
was  used  to  carry  out  this  integration  for  a  representative  set  of  values  of  9i.  The  result  was  as  expected, 
the  integral  always  evaluates  to  unity.  In  getting  this  result  it  was  necessary  to  allow  9^  to  range  over  [0  , 


Certain  commercial  equipment,  instruments,  or  materials  are  identified  in  this  paper  in  order  to  specify 
the  experimental  procedure  adequately.  Such  identification  is  not  intended  to  imply  recommendation  or 
endorsement  by  the  National  Institute  of  Standards  and  Technology  or  the  Naval  Research  Eaboratory, 
nor  is  it  intended  to  imply  that  the  materials  or  equipment  identified  are  necessarily  the  best  available  for 
the  purpose. 


7i]  rather  than  be  confined  to  the  upper  hemisphere  (the  surround)  [0, 7i/2].  The  physical  meaning  of  this 
finding  is  that  the  model,  as  outlined  above,  allows  radiation  to  be  scattered  below  the  surface.  In  fact, 
this  is  not  possible;  such  radiation  would  be  multiply  scattered  and  a  portion  would  emerge  above  the 
surface.  From  the  practical  point  of  view,  this  is  not  a  serious  problem.  For  reasonable  choices  of  (7,  the 
contribution  to  the  integral  from  the  range  [7i/2, 7i]  is  small.  In  fact,  the  light  scattered  below  the  surface 
by  the  model  can  be  thought  of  as  contributing  an  unpolarized  Lambertian  term  to  the  BRDF.  This  idea 
could  be  developed  into  a  model  for  a  depolarizing  contribution  to  the  BRDF. 

6.  Polarization 


The  discussion  to  this  point  has  been  presented  in  the  context  of  unpolarized  electromagnetic  radiation. 
However,  it  can  be  readily  generalized  to  the  case  of  polarized  radiation.  The  reason  that  the 
generalization  is  straightforward  is  that  the  microfacet  model  is  deterministic  in  nature.  Once  the  incident 
and  scattering  angles  are  specified,  the  orientation  of  the  reflecting  microfacet  is  frilly  determined.  A 
consequence  is  that  the  model  does  not  depolarize  radiation^.  If  the  incident  radiation  has  a  degree  of 
polarization,  DOP  =  1 ,  the  scattered  radiation,  at  any  scattering  angle,  will  also  have  DOP  =  1 .  This 
allows  a  Jones  matrix  approach  to  be  followed  below.  The  idea  for  using  the  Jones  matrix  approach  was 
mentioned  in  Ref.  3,  but  the  results  there  are  not  presented  in  the  useful  form  of  Eq.  (14)  below,  nor  is  the 
route  to  the  derivation  clear. 

The  generalization  of  the  material  in  Sec.  1  that  we  desire  is  to  describe  radiances  by  four-component 
Stokes  vectors  and  to  find  a  4  x  4  matrix  generalization  of  the  function  f.  This  suggests,  correctly,  that / 
should  have  the  properties  of  a  Mueller  matrix. 

For  the  purposes  of  this  section  $[,  and  are  considered  to  be  fixed.  This  also  fixes  the  orientation 
of  the  reflecting  microfacet  and  the  angles  0  and  ji  through  Eqs.  (12).  The  derivation  is  similar  to  that  of 
the  scalar  facet  model,  except  that  the  vector  nature  of  the  reflection  coefficient  p  must  be  taken  into 
account. 

The  Jones  matrix  involves  the  electric  field  vectors.  It  is  necessary  to  introduce  a  series  of  four  2- 
dimensional  (x,y)  coordinate  systems  to  track  the  electric  field  though  the  scattering  process.  The  axes  of 
these  coordinate  systems  are  the  familiar  s  and  p  of  Fresnel  reflection  theory.  The  designation  s  means 
perpendicular  to  a  reference  plane,  and  the  designation  p  means  parallel  to  the  reference  plane,  both  are 
perpendicular  to  a  propagation  direction.  For  the  first  coordinate  system,  the  reference  plane  is  defined  by 
the  incident  direction  and  the  surface  normal  z.  For  the  second  coordinate  system  the  reference  plane  is 
defined  by  the  microfacet  normal  n  and  the  incident  direction.  The  second  coordinate  system  is  rotated 
away  from  the  first  (about  the  incident  direction)  by  an  angle  rji.  It  will  be  seen  below  that  rji  is  a  function 
of  dr,  and  For  the  third  coordinate  system,  the  reference  plane  is  defined  by  the  microfacet 
normal  n  and  the  reflected  direction.  This  reference  plane  is  the  same  as  for  the  second  coordinate  system 
but  the  propagation  direction  (needed  to  specify  p)  is  now  the  reflected  direction.  For  the  fourth 
coordinate  system,  the  reference  plane  is  defined  by  the  reflected  direction  and  the  surface  normal  z.  The 
fourth  coordinate  system  is  rotated  away  from  the  third  (about  the  reflected  direction)  by  an  angle  —rjr-  In 
general,  the  angle  Pr  is  not  equal  to  rji,  but  is  also  determined  uniquely  by  0,  9^,  and 

The  Fresnel  amplitudes  for  ss  (ass)  and  pp  (app)  reflection  relate  the  reflected  electric  field,  as  described  in 
the  third  coordinate  system,  to  the  incident  electric  field,  as  described  in  the  second  coordinate  system. 
Schematically  this  can  be  represented  as  Er  =  a  Ei  where  a  is  a  complex  diagonal  2x2  matrix.  What  is 
needed  is  the  Jones  matrix,  a  2  x  2  matrix  relating  the  reflected  electric  field,  as  described  in  the  fourth 
coordinate  system,  to  the  incident  electric  field,  as  described  in  the  first  coordinate  system.  With  the  aid 
of  the  angles  rji  and  rj^,  and  rotation  matrices,  this  can  be  given  as: 


^  cos(77j  sin(;7^)^ 

0  ^ 

^cos(rip  -sin(r}p'' 

IfJ 

v-sin(/7,)  cos{riP^ 

V  0  ^PP, 

cos(;7,)  ^ 

(T 

ss 

T  ^ 

ps 

IfJ 

T 

T 

pp  J 

(14) 


The  second  relationship  is  the  working  definition  of  the  Jones  matrix,  T.  Here  and  Upp  are  Fresnel 
amplitudes,  not  reflectivities.  In  general  they  are  complex  functions  of  P  and,  through  P,  functions  of  6i, 
and  (p-cjh.  Explicit  forms  for  the  elements  of  the  Jones  matrix  in  terms  of  Uss,  flpp,  rji  and  rj^  can  be 
easily  worked  by  carrying  out  the  matrix  multiplications  above. 

To  complete  the  analysis  it  is  necessary  to  give  an  expression  relating  the  Mueller  matrix  elements  to 
those  of  the  Jones  matrix,  and  to  provide  explicit  expressions  for  rj^  and  rj^.  Reference  9  gives  a  formal 
relation  between  the  Jones  matrix,  T,  and  the  Mueller  matrix,  M.  Perhaps  more  directly  useful  are  the 
explicit  formulas  found  in  Bohren  and  Huffman^®.  Allowing  for  differences  in  notation  these  are: 


+\tJ) 

Mo2=i(7;j;+cc  +tx,+cc) 
^03  =  iW7;7:-cc)+/(7;X-cc)] 

M,,=\[iTX+cc)  -iTX+^^)] 
M,,=^[i(TX-cc)-i(TXl-cc)] 
M2,=WssK+CC  +TXp+^^) 
M,,=\[{TX+cc)  -iTX,+cc)\ 
M22=XssT;p+^^  +Xsl+^^) 
M22  =  i  ixK  -  cc)  -  ktXp  - 

^30  =W(TssK-^^)+XsKp 

A^31  =  m^ssK -CC)-i{TXl  -cc)] 

M,,=\[i{TX-cc)  +  i{TXXcc)] 
m,,=\[{tX+cc) 


(15) 


Here  the  asterisk  denotes  the  operation  of  taking  the  complex  conjugate,  cc  stands  for  complex  conjugate 
of  the  immediately  preceding  expression,  and  i  is  the  imaginary  unit. 

The  expression  for  the  polarimetric  BRDF  is  therefore 


(16) 


where  j  and  k  each  range  from  0  to  3.  In  an  expression  such  as  Eq.  (2)  ,  i?sky  and  i^reflected  are  replaced  by 
Stokes  vectors,  and  the  scalar  BRDF  is  replaced  by  the  matrix  form  above. 

Eq.  (15)  simplifies  in  some  cases.  For  the  case  of  in-plane  scattering  (^4-i^i  =0),  rji  =  rj^  =0,  and  the  off 
diagonal  elements  of  the  Jones  matrix  are  zero.  This  leaves  only  Mqo,  Mqi,  M\o,  M\\,  Mu,  M^,  M32,  and 
M33  non-zero,  with  Mqi  =  M\o,  Moo  =  Mn,  M22  =  M33  and  M23  =  -M32.  This  is  the  pattern  often  seen  in 
literature  discussions  of  the  Mueller  matrix.  The  fact  that  Moo  =  1  for  this  case  is  indicative  of  the  fact 

that  there  is  no  depolarization  in  this  model.  A  simpler  case  occurs  when  all  elements  of  the  Jones  matrix 
are  real,  which  is  the  situation  in  Fresnel  theory  when  the  index  of  refraction  is  real.  In  this  case  M03  = 
M30  =  Mi3  =  M31  =  M23  =  M32  =  0,  and  no  circular  polarization  can  be  generated  unless  the  incident 
radiation  has  a  circular  component. 

The  remaining  task  is  to  give  expressions  for  the  angles  rji  and  in  terms  of  di,  9^,  and  For  this 
purpose  it  is  useful  to  introduce  four  unit  vectors".  These  are  Si,  Zi ,  Sr,  and  2^.  As  above,  the  subscript  i 
refers  to  the  incident  ray  and  the  subscript  r  to  the  reflected  ray.  The  Si,  unit  vector  is  perpendicular  to 
the  incident  direction  and  to  the  z  direction.  It  is  parallel  to  the  electric  field  for  radiation  polarized 
perpendicular  to  the  reference  plane  in  the  first  of  the  four  coordinate  systems  defined  above.  The  Zj  unit 
vector  is  perpendicular  to  the  incident  direction  and  to  the  n  direction.  It  is  parallel  to  the  electric  field  for 
radiation  polarized  perpendicular  to  the  reference  plane  in  the  second  of  the  four  coordinate  systems. 

The  Zr  unit  vector  is  perpendicular  to  the  reflected  direction  and  to  the  n  direction.  It  is  parallel  to  the 
electric  field  for  radiation  polarized  perpendicular  to  the  reference  plane  in  the  third  coordinate  system. 
Finally,  the  Sr  unit  vector  is  perpendicular  to  the  reflected  direction  and  to  the  z  direction.  It  is  parallel  to 
the  electric  field  for  radiation  polarized  perpendicular  to  the  reference  plane  in  the  fourth  coordinate 
system.  These  four  unit  vectors  can  be  related  to  the  unit  vectors  introduced  earlier  by  the  cross  product 
relations: 


S .  =  (r.  X  z)  / 1  r.  X  z  I  =  (r.  X  z)  /  sin(0. ) 

2,.  =  (r.  X  n)  /  I  r.  X  n  I  =  (r.  x  n)  /  sin(y5) 

=  -(r^  X  z)  / 1  X  z  I  =  -  (r^  X  z)  /  sin(6>^ ) 
X  n)  / 1  X  n  I  =  -  (r^  X  n)  /  sin(y5) 
n  =  (r,+*',)/(2cos(/?)) 


The  expression  for  n  is  repeated  from  Eq.  (7)  above.  From  these  definitions  it  can  be  seen  that 
cos(;7. )  =  S,  ■  S .  =  (r.  X  z)  ■  (r.  x  n)  /(sin(6>. )  sin(^))  = 

cos(;7.)  =  S,  ■  2,  =  (r,  x  z)  ■  (r,  xn)  /(sin(0 J  sin(^))  = 

sin(6>Jsin(y») 


A  standard  vector  identity  was  used  to  eliminate  the  cross  products.  The  expression  for  n  given  in  Eq.  (7) 
can  be  substituted  into  these  expression  to  reduce  this  expression  to  trigonometric  functions  of  9^,  9^,  and 
P.  The  results  can  be  written  as: 


cos(r}. )  =  [(cos(0, )  +  cos(0^ ))  /(2  cos(^))  -  cos(0. )  cos(/?)]  /(sin(0. )  sin(/?)) 

cos(;7^ )  =  [(cos(6>, )  +  cos(6>^ ))  /(2  cos(y5))  -  cos(6>, )  cos(y5)]  /(sin(6>^ )  sin(y5))  ( 1 9) 

cos(2/?)  =  cos(0, )  cos(0^ )  +  sin(0,. )  sin(0^ )  cos(^^  -  ) 

The  equation  for  cos(2y0)  follows  from  the  spherical  addition  theorem.  These  three  equations  constitute 
the  expression  for  the  angles  rj,  and  in  terms  of  0i,  9^,  and  It  can  be  easily  checked  that  ?7i  =  ?7r  = 
0  when  and  (j),  =  (jf^. 


7.  Visualization  of  the  BRDF 


The  analytic  formulas  above  are  too  complex  to  allow  the  character  of  the  BRDF  to  be  directly 
apprehended.  Contour  plots  in  the  6„  plane  (for  fixed  incident  direction)  are  a  useful  way  to  visualize 
the  structure  of  the  BRDF.  This  section  presents  examples  of  this  type  of  plot.  The  values  of  the 
parameters  have  been  chosen  not  to  match  data  but  to  allow  certain  key  features  to  be  easily  seen.  As  a 
further  aid  to  visualization,  the  quantity  plotted  is  the  product  of  the  BRDF  component  and  cos(0r)- 
Without  multiplication  by  this  factor  of  cosine,  the  in-plane  BRDF  diverges  as  9^  approaches  90  degrees. 
This  divergence  is  real  in  the  sense  that  it  is  seen  in  measured  data  where  the  effective  value  of  (7  is  less 
than  about  1.5.  The  divergence  is  confined  to  near  in-plane  angles  and  does  not  present  problems  in 
practical  use  of  the  BRDF  because  of  the  cosine  factors  in  integration  formulas  such  as  Eq.  (1)  and  (2). 

The  parameter  values  chosen  for  Figs.  1-4  are  (7=  0.15,  and  index  of  refraction  n  =  1.57,  with  an 
incidence  angle  of  60  degrees.  Since  the  value  of  n  is  real,  the  symmetries  discussed  in  the  previous 
section  apply  and  Brewster  angle  effects  are  significant.  Contour  plots  for  the  00,  02,  22  and  33  matrix 
elements  are  shown  in  Figs.  1  through  4,  respectively.  In-plane  scattering  corresponds  to  ^  =  0  degrees  in 
these  plots.  It  can  be  noted  the  00  and  1 1  components  peak  near  the  specular  scattering  direction  as 
expected.  However,  their  behavior  out-of-plane  is  significantly  different.  The  general  behavior  of  the  0 1 
component  (not  shown)  is  similar  to  that  of  the  00  component.  The  22  component  shows  dramatic  effects 
of  passing  through  the  Brewster  angle.  The  track  of  the  Brewster  angle  position  can  be  seen  in  the  9^,  (pr 
plane.  The  02  component  is  zero  for  in-plane  scattering  as  expected  from  symmetry.  The  contour  plot 
shows  a  peak  approximately  10  degrees  removed  from  the  in-plane  condition.  The  12  component  (not 
shown)  exhibits  similar  behavior.  The  complex  behavior  of  these  off-diagonal  components  suggests  the 
attempts  to  fit  to  simple  functional  forms  will  encounter  difficulties. 
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8.  Comparison  with  measurements 


The  ideal  data  set  for  comparison  purposes  would  feature  full  polarimetric  measurements,  both  in-plane 
and  out-of-plane,  at  several  wavelengths.  Such  a  full  data  set  does  not  appear  to  exist.  We  have  obtained 
data  sets  measured  on  metal  coupons  painted  with  CARC  Green  383  paint  from  Surface  Optics 
Corporation  (SOC)'^  and  from  the  National  Institute  of  Standards  and  Technology  (NIST).  The  SOC  data 
set  is  limited  to  in-plane  data,  but  is  polarimetric  and  contains  measurements  at  three  representative 
wavelengths  (0.50  Rm,  4.60  Rm  and  9.925  |J,m).  The  NIST  data,  while  confined  to  the  visible  (0.633  Rm), 
features  (partial)  polarimetric  measurements  for  both  in-plane  and  out-of-plane  geometries.  This  section 
presents  a  preliminary  investigation  of  the  fit  between  the  microfacet  model  and  the  measured  data. 

The  model  described  above  has  just  three  parameters,  the  slope  distribution  width  parameter  a,  and  the 
real  and  imaginary  parts  of  the  index  of  refraction  n.  It  is  a  rather  strenuous  exercise  to  fit  a  tensor 
function  over  three  angles  with  just  three  parameters.  As  will  be  seen  below,  the  model  predicts  major 
trends  correctly  but  is  not  a  close  numerical  fit.  For  this  preliminary  investigation,  the  parameters  are 
used  freely  as  fit  parameters  with  only  general  reasonability  used  as  a  constraint. 

In  the  case  of  the  visible  (0.633  Rm)  measurements,  the  parameter  space  reduces  to  two  because  it  is  not 
expected  that  the  imaginary  part  of  the  index  of  refraction  would  play  a  major  role.  A  small,  but  not 
negligible,  amount  of  depolarized  scattering  is  evident  on  examining  the  NIST  data  set.  As  discussed 
above,  this  contribution  to  the  BRDF  is  not  modeled  by  the  microfacet  model.  In  the  NIST 
measurements^,  the  sample  was  illuminated  with  s-  or  p-polarized  radiation,  and  the  BRDF  was  recorded. 
A  first  order  empirical  allowance  for  the  depolarized  scattering  can  be  made  by  introducing  a  diffuse 
contribution  additive  to  the  model  BRDF.  In  the  figures  below  this  is  taken  as  a  constant  additive 
contribution  of  0.025  sr“^  to  the  BRDF  for  the  case  of  j3-polarized  incident  radiation,  and  a  contribution  of 
0.04  sr“’  for  s-polarized  incident  radiation. 
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Figure  6 


^The  uncertainties  in  the  data  are  dominated  by  statistical  sources,  and  can  be  estimated  by  the  point-to- 
point  spread  in  the  data.  Other  sources  of  uncertainty  are  expected  to  be  below  3  %  of  any  given  value. 


Figure  5  shows  the  in-plane  BRDF  for  an  incident  angle  of  65  degrees,  shown  as  a  function  of  polar 
scattering  angle.  The  model  parameters  are  as  shown  in  the  text  box  in  the  figure.  The  theoretical  curves 
are  (^o+/ol)  +  0.04  for  the  case  of  s-polarization  and  (fm-fm)  +  0.025  for  the  case  of p-polarization.  As 
can  be  seen  from  the  curves,  the  qualitative  features  of  the  data  are  captured  by  the  model.  Light  which  is 
s-polarized  is  always  more  strongly  scattered  than  light  which  is  p-polarized,  and  there  is  no  discernable 
peak  at  the  specular  scattering  angle  of  65  degrees.  The  BRDF  diverges  as  the  polar  scattering  angle 
approaches  90  degrees.  Flowever,  there  is  no  combination  of  the  two  parameters  (n  and  o)  which  match 
the  BRDF  over  the  whole  angular  range.  For  the  parameter  choice  shown,  the  divergence  is  weaker  in  the 
data  than  predicted  by  the  model.  Further  the  value  n=  1.25  is  lower  than  might  have  been  expected  on 
physical  grounds. 

Figure  6  shows  the  out-of-plane  BRDF  for  p-polarized  radiation,  shown  as  a  function  of  azimuthal  angle 
for  fixed  polar  angles  (both  incident  and  scattering)  of  65  degrees.  As  with  the  case  of  the  in-plane  data, 
the  qualitative  features  are  captured  but  quantitative  agreement  is  not  obtained.  The  wide  width  seen  in 
the  data  indicates  that  the  functional  form  of  the  actual  slope  distribution  falls  off  more  gently  than  the 
exponential  form  used  in  the  model. 

The  NIST  data  also  allows  the  degree  of  polarization  (DOP)  to  be  studied.  For  the  out-of-plane  data  of 
Fig.  6,  the  DOP  decreases  as  the  azimuthal  angle  increases  from  about  0.7  to  about  0.4.  This  is  roughly 
consistent  with  the  model  predictions  assuming  that  the  0.025  additive  contribution  is  approximately 
correct. 

A  different  test  of  the  polarization  performance  of  the  microfacet  model  can  be  made  by  an  examination 
of  the  orientation  of  the  principal  direction  of  linear  polarization,  as  the  azimuthal  scattering  angle  is 
varied.  For  the  data  of  Fig.7  the  incident  radiation  is  p-polarized  and  both  the  incident  and  scattering 
polar  angles  are  set  at  45  degrees.  The  scattering  azimuth  is  varied,  and  the  orientation  of  the  principal 
direction  is  plotted.  A  value  of  90  degrees  for  this  angle  means  that  the  scattered  linearly  polarized  light 
is  j3-polarized,  while  a  value  of  0  degrees  corresponds  to  s-polarization.  As  can  be  seen  from  the  Fig.  7, 
the  scattered  light  is  p-polarized  for  in-plane  scattering  (as  expected  from  symmetry)  and  tends  to  be  s- 
polarized  as  the  scattering  azimuth  increases  away  from  the  in-plane  configuration.  A  very  good  fit  of 
the  microfacet  model  to  the  data  was  possible  with  a  choice  ofn  =  1.5. 


Figure  7 

Figures  8  and  9  present  an  in-plane  BRDF  measured  at  SOC  using  9.925  gm  radiation.  In  Fig.  8,  the /oo 
component,  for  incident  angle  equal  to  60  degrees,  is  plotted  as  a  function  of  polar  scattering  angle.  The 
model  parameters  are  shown  in  the  figure.  The /lo  component  is  shown  in  Fig.  9.  For  the  9.925  gm 
radiation  case  it  is  necessary  to  include  an  imaginary  part  of  the  index  of  refraction  to  get  qualitative 


agreement  between  measurements  and  the  model.  However,  the  depolarized  diffuse  contribution  is 
relatively  smaller  and  can  be  neglected  for  purposes  of  this  preliminary  analysis. 
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Figure  8  Figure  9 

As  can  be  seen  in  the  figures,  qualitative  features  of  the  data  are  captured  by  the  model.  There  is  a  peak 
near  the  specular  scattering  angle  (60  degrees)  which  is  shifted  to  larger  angle  by  the  (weak)  divergence 
as  90  degrees  is  approached.  The  numerical  values  are  in  the  correct  range  until  past  the  specular  angle. 
The  modulus  of  the  complex  index  of  refraction  chosen  to  obtain  the  fit  is  lower  than  would  be  expected 
on  physical  grounds. 

The  fit  value  of  the  model  width  parameter,  a,  is  only  0.1  for  the  9.925  |J,m  case  as  compared  with  0.3  for 
the  visible  case.  In  the  model  calculations,  this  is  why  the  9.925  |J,m  BRDF  shows  a  peak  at  the  specular 
scattering  angle  while  the  visible  BRDF  does  not.  It  is  reasonable  to  ascribe  the  difference  in  behavior  at 
the  two  wavelength  regions  to  diffraction  effects  not  directly  modeled  in  the  microfacet  model.  A  non¬ 
rigorous  argument  can  be  advanced  to  motivate  the  use  of  a  smaller  value  of  cr  for  the  9.925  |J,m  case 
than  for  the  0.633  |J,m  case.  The  value  of  is  the  value  of  the  surface  slope  variance.  The  variance  can 
be  viewed  as  the  integral,  over  wavenumber,  of  the  power  spectral  density  of  surface  slope.  To  perform 
the  integral,  an  upper  wavenumber  cut-off  must  be  selected.  The  higher  the  cut  off,  the  larger  the 
resulting  variance.  It  is  reasonable  to  argue  that  diffraction  effects  should  set  the  cut  -off  at  roughly  liilX 
where  X  is  the  wavelength  of  the  light  being  scattered.  This  argument  suggests  that  the  cut  off  for  0.633 
pm  light  is  some  20  times  higher  than  for  the  9.925  pm  light 


9.  Summary 


Closed  form  expressions  for  the  full  polarimetric  form  of  the  BRDF  has  been  derived  for  the  microfacet 
model.  Numerical  study  of  this  result  reveals  complex  structure  for  some  of  the  individual  components  of 
the  BRDF,  especially  the  off-diagonal  elements.  Comparisons  with  measured  data  indicate  that  the  three 
parameter  model  captures  important  qualitative  features  of  the  data  but  falls  well  short  of  numerical 
agreement. 
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